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Abstract. The exclusive uj electroproduction off the proton was studied in a large kinematical domain 
above the nucleon resonance region and for the highest possible photon virtuality (Q 2 ) with the 5.75 GeV 
beam at CEBAF and the CLAS spectrometer. Cross sections were measured up to large values of the 
four-momentum transfer (— t < 2.7 GeV 2 ) to the proton. The contributions of the interference terms <jtt 
and otl to the cross sections, as well as an analysis of the u spin density matrix, indicate that helicity is 
not conserved in this process. The i-channel ty° exchange, or more generally the exchange of the associated 
Regge trajectory, seems to dominate the reaction 7*p — > cop, even for Q 2 as large as 5 GeV 2 . Contributions 
of handbag diagrams, related to Generalized Parton Distributions in the nucleon, are therefore difficult to 
extract for this process. Remarkably, the high-£ behaviour of the cross sections is nearly Q 2 -independent, 
which may be interpreted as a coupling of the photon to a point-like object in this kinematical limit. 

PACS. 13.60.Le Production of mesons by photons and leptons - 12.40.Nn Regge theory - 12.38.Bx 
Perturbative calculations 
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1 Introduction 

The exclusive electroproduction of vector mesons is a pow- 
erful tool, on one hand to understand the hadronic prop- 
erties of the virtual photon (7*) which is exchanged be- 
tween the electron and the target nucleon p^, and on the 
other hand to probe the quark-gluon content of the pro- 
ton (p) I21CI1E]- At moderate energies in the -y*p system, 
but large virtuality of the photon, quark-exchange mech- 
anisms become significant in the vector meson production 
reactions j*p — > pp/u>, thus shedding light on the quark 
structure of the nucleon. 

The interaction of a real photon with nucleons is dom- 
inated by its hadronic component. The exchange in the 
i-channel of a few Regge trajectories permits a descrip- 
tion of the energy dependence as well as the forward an- 
gular distribution of many, if not all, real-photon-induced 
reactions (see e.g. Ref. For instance, this approach re- 
produces the photoproduction of vector mesons from the 
CEBAF energy range to the HERA range (a few to 200 
GeV) 0. The exchange of the Pomeron (or its realiza- 
tion into two gluons) dominates at high energies, while 
the exchange of meson Regge trajectories (it, <j, f 2 ) takes 
over at low energies. At jp energies of a few GeV, lu pho- 
toproduction off a proton is dominated by tt° exchange 



in the i-channel (fig. [TJ. The use of a saturating Regge 
trajectory 0J is very successful in describing recent pho- 
toproduction data [7] at large angles (large momentum 
transfer t). This is a simple and economical way to pa- 
rameterize hard scattering mechanisms. Extending these 
measurements to the virtual photon sector opens the way 
to tune the hadronic component of the exchanged pho- 
ton, to explore to what extent 7r° exchange survives, and 
to observe hard scattering mechanisms with the help of a 
second hard scale, the virtuality Q 2 of the photon. 

The study of such reactions in the Bjorken regime 1 
holds promise, through perturbative QCD, to access the 
so-called Generalized Parton Distributions (GPD) of the 
nucleon [8||9]. These structure functions are a generaliza- 
tion of the parton distributions measured in the deep in- 
elastic scattering experiments and their first moment links 
them to the elastic form factors of the nucleon. Their sec- 
ond moment gives access to the sum of the quark spin and 
the quark orbital angular momentum in the nucleon 
The process under study may be represented by the so- 
called handbag diagram (fig.GJ. Its amplitude factorizes [TH] 



1 Q 2 and v large and xb finite, where —Q 2 and v are the 
squared mass and the laboratory-frame energy of the virtual 
photon, while xb = Q 2 /2M p u is the usual Bjorken variable. 




Fig. 1. Schematic representations of the t-channel exchange 
(left) and of the handbag diagram (right) for exclusive vector 
meson electroproduction. 



into a "hard" process where the virtual photon is ab- 
sorbed by a quark and a "soft" one containing the new 
information on the nucleon, the GPD (which are func- 
tions of x and x' , the momentum fraction carried by the 
quark in the initial and final states, and of t, the squared 
four-momentum transfer between the initial and final pro- 
tons). The factorization applies only to the transition, 
at small values of —t, between longitudinal photons (L) 
and helicity-0 mesons, which is dominant in the Bjorken 
regime. Because of the necessary gluon exchange to pro- 
duce the meson in the hard process (see fig. 0) , the domi- 
nance of the handbag contribution is expected to be reached 
at a higher Q 2 for meson production than for photon pro- 
duction (DVCS). Nevertheless, recent results on deeply 
virtual p production show a qualitative agreement with 
calculations based on the handbag diagram [111112) . Vector 
meson production is an important complement to DVCS, 
since it singles out the quark helicity independent GPD H 
and E which enter Ji's sum rule t 8 and allows, in princi- 
ple, for a flavor decomposition of these distributions (see 
e.g. Ref. EH). 

Apart from early, low statistics, muon production ex- 
periments at SLAC |14U15| . the leptoproduction of u> mesons 
was measured at DESY QH, for 0.3 < Q 2 < 1.4 GeV 2 , 
W < 2.8 GeV (x B < 0.3), and then at Cornell QH, in a 
wider kinematical range (0.7 < Q 2 < 3 GeV 2 , W < 3.7 
GeV) but with larger integration bins. These two experi- 
ments yielded cross sections differing by a factor of about 2 
wherever they overlap (around Q 2 ~ 1 GeV 2 ). The DESY 
experiment also provided the only analysis so far, in elec- 
troproduction, of the lo spin density matrix elements, av- 
eraged over the whole kinematical range. This analysis in- 
dicated that, in contrast with p electroproduction, there is 
little increase in the ratio R of longitudinal to transverse 
cross sections (ul/o't) when going from photoproduction 
to low Q 2 electroproduction. More recently, lo electropro- 
duction was measured at ZEUS ^H], a * high Q 2 and very 
low x b , in a kinematical regime more sensitive to purely 
diffractive phenomena and to gluons in the nucleon. Fi- 
nally, there is also unpublished data from HERMES |19| . 

The main goal of the present experiment was to reach 
the highest achievable Q 2 values in exclusive meson elec- 
troproduction in the valence quark region. In the specific 
case of the lo production, it is to test which of the two de- 
scriptions — with hadronic or quark degrees of freedom, 
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Fig. 2. (Color online) Kinematical range covered by this and 
previous |16I17I18I19| u> electroproduction experiments. The 
lines are indicative of the total coverage in Q 2 and xg of pre- 
vious experiments. 



more specifically i-channel Regge trajectory exchange or 
handbag diagram — applies in the considered kinematical 
domain (see fig. |5J. For this purpose, the reduced cross 
sections cr 7 * p ^ wp were measured in fine bins in Q 2 and 
xb, as well as their distribution in t and <fi (defined be- 
low). In addition, parameters related to the to spin density 
matrix were extracted from the analysis of the angular dis- 
tribution of the lo decay products. If the vector meson is 
produced with the same helicity as the virtual photon, 
s-channel helicity conservation (SCHC) is said to hold. 
From our results, the relevance of SCHC and of natural 
parity exchange in the t-channel was explored in a model- 
independent way. These properties have been established 
empirically in the case of photo- and electroproduction of 
the p meson (see e.g. Ref. HO]), but may not be a general 
feature of all vector meson production channels. 

This paper is based on the thesis work of Ref. |21) . 
where additional details on the data analysis may be found. 



2 Experimental procedure 

We measured the process ep — > epto, followed by the decay 
lo — > 7T + 7r~7r°. The scattered electron and the recoil pro- 
ton were detected, together with at least one charged pion 
from the lo decay. At a given beam energy E, this process 
is described by ten independent kinematical variables. In 
the absence of polarization in the ep initial state, the ob- 
servables are independent of the electron azimuthal angle 
in the laboratory. Q 2 and xb are chosen to describe the 
7*p initial state. The scattered electron energy E 1 and, 
for ease of comparison with other data, the 7*p center- 
of-mass energy W will be used as well, t is the squared 
four- momentum transfer from the 7* to the lo, and <p the 
angle between the electron (ee'7*) and hadronic {"f*Lop) 
planes. Since t is negative and has a kinematical upper 
bound t (Q 2 , Xb) corresponding to lo production in the di- 
rection of the 7*, the variable t' = to — t will also be used. 
The lo decay is described in the so-called helicity frame, 
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Fig. 3. (Color online) Reference frames and relevant variables 

used for the description of the reaction ep — > epuj, followed by 
+ - o 

U! — > 7T 7T 7T . 



where the w is at rest and the z-axis is given by the lu 
direction in the "f*p center-of-mass system. In this helicity 
frame, the normal to the decay plane is characterized by 
the angles On and tp N (fig.|3J|. Finally the distribution of 
the three pions within the decay plane is described by two 
angles and a relative momentum. This latter distribution 
is known from the spin and parity of the u> meson |22| 
and is independent of the 7*_p — > top reaction mechanism. 
The purpose of the present study is to characterize as 
completely as possible the distributions of cross sections 
according to the six variables Q 2 , xb, t, 4>, cos On and (fiN- 



2.1 The experiment 

The experiment was performed at the Thomas Jefferson 
National Accelerator Facility (JLab). The CEBAF 5.754 
GeV electron beam was directed at a 5-cm long liquid- 
hydrogen target. The average beam intensity was 7 nA, 
resulting in an effective integrated luminosity of 28.5 fb _1 
for the data taking period (October 2001 to January 2002). 
The target was positioned at the center of the CLAS spec- 
trometer. This spectrometer uses a toroidal magnetic field 
generated by six superconducting coils for the determina- 
tion of particle momenta. The field integral varied approx- 
imately from 2.2 to 0.5 Tm, in average over charges and 
momenta of different particles, for scattered angles be- 
tween 14° and 90°. All the spectrometer components are 
arranged in six identical sectors. Charged particle trajec- 
tories were detected in three successive packages of drift 
chambers (DC), the first one before the region of magnetic 
field (Rl), the second one inside this region (R2), and the 
third one after (R3). Threshold Cerenkov counters (CC) 
were used to discriminate pions from electrons. Scintilla- 
tors (SC) allowed for a precise determination of the par- 
ticle time-of- flight. Finally, a segmented electromagnetic 
calorimeter (EC) provided a measure of the electron en- 
ergy. This geometry and the event topology are illustrated 




Fig. 4. (Color online) Schematic view of the CLAS spectrome- 
ter components (see text for description) and of typical particle 
tracks, viewed in projection. The torus coils are not shown. 



in fig. 01 A detailed description of the CLAS spectrometer 
and of its performance is given in Ref. |23j . 

The data acquisition was triggered by a coincidence 
CC-EC corresponding to a minimal scattered electron en- 
ergy of about 0.575 GeV. The trigger rate was 1.5 kHz, 
with a data acquisition dead time of 6%. A total of 1.25 x 
10 9 events was recorded. 



2.2 Particle identification 

After calibration of all spectrometer subsystems, tracks 
were reconstructed from the DC information. The identi- 
fication of particles associated with each track proceeded 
differently for electrons and hadrons. 

Electrons were identified from the correlation between 
momentum (from DC) and energy (from EC). In addi- 
tion pions were rejected from the electron sample by a 
cut in the CC amplitude and imposing a condition on 
the energy sharing between EC components compatible 
with the depth profile of an electromagnetic shower. Ge- 
ometrical fiducial cuts ensured that the track was inside 
a high efficiency region for both CC and EC. The effi- 
ciencies of the electron identification cuts (rjee and tjec) 
depended on the electron momentum and angle (or on 
Q 2 and xb)- Vec was calculated from data samples using 
very selective CC cuts in order to unambiguously select 
electrons. r\cc was extracted from an extrapolation of the 
CC amplitude Poisson distribution into the low amplitude 
region. These efficiencies varied respectively between 0.92 
and 0.99 (CC), and 0.86 to 0.96 (EC). At low electron en- 
ergies, a small contamination of pions remained, which did 
not however satisfy the lo selection criteria to be described 
below. 

The relation between momentum (from DC) and ve- 
locity (from path length in DC and time-of-flight in SC) 
allowed for a clean identification of protons (p) and pi- 
ons (tt + and 7r _ ). However, for momenta larger than 2 
GeV/c, ambiguities arose between p and ir + identification, 
which led us to the discarding of events corresponding to 
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Fig. 5. (Color online) Identification of the uj channel in the 
case of the detection of one charged pion: M x [epn + X] vs 
Mx[epX]. Events corresponding to p (horizontal locus) and 
other two-pion production channels are clearly separated from 
those corresponding to uj production (vertical locus). 




Fig. 6. (Color online) Missing mass Mx[epX] distributions for 
the ep — ► epir + X event configuration, for two (Q 2 ,xb) bins, 
after selection cuts and event weighting discussed in the text. 
Left : 2.2 GeV 2 < Q 2 < 2.5 GeV 2 and 0.34 < x B < 0.40. 
Right : 3.1 GeV 2 < Q 2 < 3.6 GeV 2 and 0.52 < x B < 0.58 (at 
the edge of kinematical acceptance) . The two lines indicate the 
subtracted background (green) and the fitted distribution (a; 
peak + background in red). 



t < —2.7 GeV 2 . Fiducial cuts were applied to hadrons 
as well. The efficiency for the hadron selection cuts was 
accounted for in the acceptance calculation described in 
sect. El 

2.3 Event selection and background subtraction 

Two configurations of events were studied, with one or two 
detected charged pions: ep — > epn + X and ep — > epir + ir~ X. 
The former benefits from a larger acceptance and is ade- 
quate to determine cross sections, while the latter is nec- 
essary to measure in addition the distribution of the uj 
decay plane orientation and deduce from it the uj spin 
density matrix. The final selection of events included cuts 
in W and E': W > 1.8 GeV to eliminate the threshold 
region sensitive to resonance production [211 an d E' > 0.8 
GeV to minimize radiative corrections and residual pion 
contamination in the electron tracks. The first configu- 
ration was selected requiring a missing mass M x larger 
than 0.316 GeV to eliminate two pion production chan- 
nels (M| > 0.1 GeV 2 on the vertical axis of fig. [SJ . This 
cut was chosen slightly above the two pion mass in or- 
der to minimize background. The corresponding losses in 
ep — > epuj events were very small and accounted for in the 
acceptance calculation to be discussed below. Events cor- 
responding to the uj production appear as a clear peak in 
the ep — > epX missing mass spectrum (fig. |H} . The width 
of this peak (er ~ 16 MeV) is mostly due to the experi- 
mental resolution. 

After proper weighting of each event with the accep- 
tance calculated as indicated in sect. 12.41 a background 
subtraction was performed for each of 34 bins (Q 2 ,Xb) 
and, for differential cross sections, for each bin in t or (f>. 
The background was determined by a fit to the acceptance- 
weighted distributions with a second-order polynomial and 
a peak shape as modeled by simulations (a skewed gaus- 
sian shape taking into account the experimental resolution 
and radiative tail) . At the smallest values of W, the fitted 
background shape was modified to account for kinemati- 
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Mx [e pX] (GeV) 

Fig. 7. (Color online) Identification of the uj channel in the 
case of the detection of two charged pions: M x [epn + ir~X] vs 
M x [epX] for M x [epTT ± X] > 0.1 GeV 2 . The spot at ( 0.78, 
0.02) corresponds to the uj. 



cal acceptance cuts. The acceptance-weighted numbers of 
ep — > epuj events were computed using the sum of weighted 
counts in the M x [epX] distributions for .72 < M x < .85 
GeV, diminished by the fitted background integral in the 
same interval. 

Likewise, events from the second configuration (ep — ► 
epTT + ir~ X) were selected with cuts in missing masses: 
M x [epir+X] and M x [ejyK~X] > 0.316 GeV, 0.01 GeV 2 < 
M x [epir+ir-X] < 0.045 GeV 2 (fig. 01. The resulting 
M x [epX) spectrum after these cuts is illustrated in fig. [8] 
The background subtraction in the spectrum of weighted 
events proceeded in the same way for each of 64 (Q 2 , xb, 
cos 8 n) bins or for each of 64 (Q 2 , xb, <Pn) bms in order 
to analyze the w decay distribution (see sect.QJ. 

2.4 Acceptance calculation 

The tracks reconstruction and the event selection were 
simulated using a GEANT-based Monte Carlo (MC) sim- 
ulation of the CLAS spectrometer. We used an event gen- 
erator tuned to reproduce photoproduction and low Q 2 
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Fig. 8. Unweighted Mx[epX] spectrum for all ep — * epir + -K~X 
events after selection cuts discussed in the text. 
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Fig. 9. 4> dependence of calculated acceptance, integrated over 
other kinematical variables, for the one and two detected pion 
configurations: ep — > epix + X (left), ep — > epix + -K~X (right). 



data in the resonance region and extrapolated into our 
kinematical domain |25| . The acceptance was defined in 
each elementary bin in all relevant variables as the ratio 
of accepted to generated MC events. At the limit of small 
six-dimensional bins, it is independent of the model used 
to generate the MC events. The MC simulation included 
a tuning of the DC and SC time resolutions to reproduce 
the observed widths of the hadron particle identification 
spectra and of the missing mass spectra, so that the effi- 
ciency of the corresponding cuts described above could be 
correctly determined. 

For the extraction of cross sections from the ep — > 
epir + X configuration, acceptance calculations were per- 
formed in 1837 four-dimensional bins (Q 2 , xb, t and (f>) 
with two different assumptions about the event distribu- 
tion in cos 6 7v and tp N . The two different MC calcula- 
tions were used for an estimate of the corresponding sys- 
tematic uncertainties (see sect. 13. If) . For the analysis of 
the decay plane distribution W(cos <pn, <fi) from the 
ep — > epTr + Tr~ X configuration, the acceptance calculation 
was performed in 3575 six-dimensional bins (Q 2 , xb, t, (j), 
cos On and ipx)- The binning is defined in table ^and the 
numbers above correspond to kinematically allowed bins 
that have significant statistics. 

The calculated acceptances are, on average, of the or- 
der of 2% and 0.2% respectively for the two event con- 
figurations of interest. They vary smoothly for all vari- 
ables except 4>, where oscillations, due to the dead zones 
in the CLAS sectors, reproduce the physical distributions 
of events (fig. El- Each event was then weighted with the 
inverse of the corresponding acceptance. Events belong- 
ing to bins with either very large or poorly determined 
weights were discarded (for the ep — > epir + X configura- 
tion, acceptance smaller than 0.25% or associated MC sta- 
tistical uncertainty larger than 35%). The corresponding 
losses (a few percent) were quantified through the MC ef- 
ficiency t\mc by applying these cuts to MC events and 
computing the ratio of weighted accepted MC events to 
generated events. No attempt was made to calculate the 
acceptance for the non-resonant three-pion background, 
so that the background shape in fig. differs from the 
physical distribution da/dMx when Mx differs from the 
to mass. 



2.5 Radiative corrections 

Radiative corrections were calculated following Ref. 
They were dealt with in two separate steps. The MC ac- 
ceptance calculation presented above took into account ra- 
diation losses due to the emission of hard photons, through 
the application of the cut Mx[epX] < 0.85 GeV. Correc- 
tions due to soft photons, and especially the virtual pro- 
cesses arising from vacuum polarization and vertex cor- 
rection, were determined separately for each bin in (Q 2 , 
xb, <$)■ The same event generator employed for the com- 
putation of the acceptance was used, with radiative effects 
turned on and off, thus defining a corrective factor F rac i- 
The ^-dependence of F ra d is smaller than all uncertainties 
discussed in sect. 13. II and was neglected. 



3 Cross sections for 7*p — > up 

The total reduced cross sections were extracted from the 
data through : 



(Q 2 ,x B ,E) 



1 



n w (Q 2 ,x B ) 



r v (Q 2 ,x Bl E) BC mt AQ 2 ■ Ax B 

^rad 



VCC TjEC TJMC 



(i) 



The Hand convention [27J was used for the definition 
of the virtual transverse photon flux ry, which includes 
here a Jacobian in order to express the cross sections in 
the chosen kinematical variables : 



r v (Q 2 ,x B ,E) 



a Q 2 1 - xb 1 
8^ M 2 E 2 xl l-e'' 



(2) 



with the virtual photon polarization parameter being de- 
fined as : 

£ = -i , Q 2 + (-E-W ' ^ 
1 T Z AEE'-Q' 2 

In eq. QJ, n w (Q 2 , Xb) is the acceptance- weighted num- 
ber of ep — > epuj events after background subtraction. 
The branching ratio of the u> decay into three pions is 
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Table 1. Definition of binning for cross section (1) and u 
polarization (2) data. N refers to the number of bins in the 
specified range for each variable. 



Table 2. Point-to-point and normalization systematic uncer- 
tainties, for integrated and differential cross sections. 



Variable 


Ranj 


;e(l) 


N(l) 


Ranj 


;e(2) 


N(2) 


Q 2 (GeV 2 ) 


1.6 - 


3.1 


5 


1.7 - 


■ 4.1 


4 




3.1 - 


5.1 


4 


4.1 - 


■ 5.2 


1 


x B 


0.16 - 


0.64 


8 


0.18 - 


■ 0.62 


4 


-t (GeV 2 ) 


0.1 - 


1.9 


6 


0.1 - 


2.1 


4 




1.9 - 


2.7 


1 


2.1 - 


■ 2.7 


1 


(rd) 


- 


2tt 


9 


- 


2tt 


6/9/12 


cos On 








-1 


- 1 


8 


(p N (rd) 








- 


2tt 


8 



B = 0.891 [2E]- The integrated effective lumimosity jC int 
includes the data acquisition dead time correction. AQ 2 
and Axb are the corresponding bin widths; for bins not 
completely filled (because of W or E' cuts on the electron, 
or of detection acceptance), the phase space AQ 2 ■ Axb 
includes a surface correction and the Q 2 and Xb central 
values are modified accordingly. The radiative correction 
factor and the various efficiencies not included in the MC 
calculation were discussed in previous sections. 

Differential cross sections in t or <p were extracted in a 
similar manner. Cross section data and corresponding MC 
data for the acceptance calculation were binned according 
to table HI 



3.1 Systematic uncertainties 

Systematic uncertainties in the cross section measurements 
arise from the determination of the CLAS acceptance, of 
electron detection efficiencies, of the luminosity, and from 
the background subtraction. They are listed in tableland 
discussed hereafter. 

Errors in the acceptance calculation may be due to 
inhomogeneity in the detectors response, such as faulty 
channels in DC or SC, to possible deviations between 
experimental and simulated resolutions in spectra where 
cuts were applied, to the input of the event generator 
(both in cross section and in decay distribution W gen ), 
to radiative corrections and finally to the event weight- 
ing procedure. The most significant of these uncertainties 
(8%) was quantified by performing a separate complete 
MC simulation varying inputs for the parameters describ- 
ing the decay distribution W gen - 

Systematic uncertainties on the electron detection effi- 
ciencies were estimated with experimental data, by vary- 
ing the electron selection cuts or the extrapolated CC am- 
plitude distribution (see sect. 12. 2J) . 

Systematic background subtraction uncertainties were 
estimated by varying the assumed background functional 
shapes. In particular, the background curvature under the 
u> peak was varied between extreme values compatible 
with an equally good fit to the distributions in fig. H3 For 
bins corresponding to low values of W, the acceptance cut 



Source of uncertainty 



da /dt da I dcj> 



CLAS acceptance 

- inhomogeneities 

- resolutions 

l (Q 2 ,x B ,t) 



gen \ 

- W 9e „ (COS 6>jV, fN, <t>) 

- radiative corrections 

- binning 

- T}MC 

Electron detection 

- Vcc 

- r\EC 

Background subtraction 



6% 
2% 
5% 
8% 
4% 
5% 
4% 

1.5% 

2% 

7-11% 



6% 
2% 
5% 



5% 
2-7% 



6% 
2% 

2-20% 



Point-to-point 



16-18% 13-14% 7-21% 



Normalization 



9-12% 14-16% 



to the right of the Mx[epX] peak induced an additional 
uncertainty. 

Finally, overall normalization uncertainties were due to 
the knowledge of target thickness (2%) and density (1%) 
and of beam integrated charge (2%). 

Errors contributing point-to-point and to the overall 
normalization are separately added in quadrature in ta- 
ble [21 For t or cf> distributions, the same uncertainties ap- 
ply, but may contribute to the overall normalization uncer- 
tainty instead of point-to-point. For example, the shape of 
the Mx [epX] distributions depends mostly on Q 2 and xb , 
not on t and 0; the background subtraction uncertainties 
are then considered as a normalization uncertainty for the 
da / dt and da/d(f> distributions. The uncertainties on r\uc 
are largest for t and <$> bins with the smallest acceptance 
(lowest and highest t values, as well as <f> ~ 180°). 



3.2 Integrated reduced cross sections 

Results for (Q 2 , xb) are given in tableland fig. EH 

For the purpose of comparison with previous data, fig. 1111 
shows cross sections as a function of Q 2 for fixed, approx- 
imately constant, values of W. When comparing different 
data sets, note that a — <jt + £&l depends on the beam 
energy through e. However, as will be shown, it is likely 
that the difference of longitudinal contributions between 
two different beam energies (£2 — £i)cl i s much smaller 
than the total cross section a. In addition, the range of 
integration in t is different for all experiments, larger in 
this work, but most of the total cross section comes from 
small —t values. A direct comparison of the cross sections 
is then meaningful. 

There is no direct overlap between the present data 
and the DESY data ^B] , but they seem to be compatible 
with a common trend. The Cornell data |17| are roughly a 
factor 2 lower than ours. Where they overlap, the Cornell 
data are also a factor 2 lower than the DESY data. We 
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Table 3. Cross sections a = ot + eul and interference terms <jtt and gtl for the reaction 7*p — > wp, integrated over 
—2.7 GeV 2 < t < to- Slope b of da/dt for —1.5 GeV 2 < t < to- Quoted uncertainties are obtained from the addition in 
quadrature of statistical uncertainties and of point-to-point systematic uncertainties as discussed in sect. 13.11 



Xb 


Q 


W 




£ 


to 


a±Aa 




O TL i A(JTL 




b±db 






(GeV 2 ) 


(GeV) 






(GeV 2 ) 


(nb) 


(nb) 


(nb) 




(GeV 




; ) 


0.203 


1.725 


2.77 





.37 


-0.09 


536 ± 96 


60 ± 87 


2 ± 30 


2 


.44 


± 


rt 

0. 


.18 


0.250 


1.752 


2.48 





.59 


-0.15 


661 ± 118 


156 ± 61 


-35 ± 24 


1 


.93 


± 


0. 


16 


0.252 


2.042 


2.63 





.43 


-0.14 


421 ± 75 


104 ± 47 


-18 ± 18 


2 


.28 


± 


0. 


16 


0.265 


2.320 


2.70 





.32 


-0.14 


344 ± 62 


58 ± 74 


-14 ± 26 


1 


.88 


± 


0. 


17 


0.308 


1.785 


2.21 





.72 


-0.25 


1139 ± 205 


310 ± 122 


-175 ± 60 


1 


.23 


± 


0. 


.17 


0.310 


2.050 


2.33 


o 


.63 


-0.23 


551 ± 98 


121 ± 42 


-66 ± 20 


1 


.90 


± 


0. 


.16 


0.310 


2.350 


2.47 





.50 


-0.21 


395 ± 71 


103 ± 41 


-49 ± 17 


2 


.0:3 


± 


0. 


16 


0.313 


2.639 


2.58 


o 


.37 


-0.20 


287 ± 52 


111 ± 48 


-9 ± 18 


1 


.90 


± 


n 
0. 


.17 


0.327 


2.914 


2.62 


o 


.28 


-0.22 


226 ± 43 


138 ± 84 


-46 ± 27 


1 


.87 


± 


0. 


.23 


0.370 


2.050 


2.09 


o 


.74 


-0.37 


1002 ± 180 


91 ± 75 


-25 ± 39 





.97 


± 


0. 


.17 


0.370 


2.350 


2.21 


o 


.65 


-0.34 


581 ± 104 


150 ± 48 


-36 ± 23 


1 


.35 


± 


0. 


.17 


0.370 


2.650 


2.32 


o 


.55 


-0.31 


380 ± 68 


85 ± 40 


-33 ± 17 


1 


.62 


± 


0. 


17 


0.370 


2.950 


2.43 


o 


.43 


-0.30 


273 ± 49 


95 ± 42 


-44 ± 17 


1 


.93 


± 


0. 


18 


0.378 


3.295 


2.51 


o 


.31 


-0.30 


230 ± 42 


17 ± 55 


-27 ± 19 


1 


.46 


± 


0. 


.18 


0.429 


2.055 


1.90 


o 


.81 


-0.59 


2203 ± 348 


54 ± 158 


143 ± 85 


1 


.03 


± 


0. 


25 


0.430 


2.350 


2.00 


o 


.74 


-0.53 


1013 ± 182 


181 ± 82 


-102 ± 43 


o 


.78 


± 


0. 


22 


0.430 


2.650 


2.10 


o 


.67 


-0.49 


626 ± 113 


90 ± 63 


-41 ± 32 


o 


.81 


± 


0. 


22 


0.430 


2.950 


2.19 


o 


.58 


-0.46 


427 ± 78 


56 ± 48 


-24 ± 23 


o 


.92 


± 


0. 


.23 


0.430 


3.350 


2.31 


o 


.45 


-0.43 


265 ± 48 


105 ± 37 


-12 ± 15 


1 


.34 


± 


0. 


.23 


0.436 


3.807 


2.41 


o 


.30 


-0.42 


191 ± 35 


128 ± 60 


-54 ± 21 


1 


,4( 


1 

± 


0. 


24 


0.481 


2.371 


1.85 


o 


.79 


-0.79 


1660 ± 265 


-34 ± 173 


241 ± 97 


o 


.92 


± 


0. 


50 


0.490 


2.651 


1.91 





.74 


-0.75 


1113 ± 177 


291 ± 109 


15 ± 62 





.68 


± 


0. 


.37 


0.490 


2.950 


1.99 





.68 


-0.69 


644 ±116 


174 ± 69 


-154 ± 35 





.23 


± 


0. 


.34 


0.490 


3.350 


2.09 





.58 


-0.64 


397 ± 72 


84 ± 49 


-49 ± 24 


1 


.39 


± 





.25 


0.490 


3.850 


2.21 





.43 


-0.60 


272 ± 50 


72 ± 46 


5 ± 20 


1 


.41 


± 





.26 


0.494 


4.307 


2.30 





.29 


-0.58 


187 ± 37 


169 ± 79 


-14 ± 27 





.78 


± 


0. 


.32 


0.538 


2.968 


1.85 





.73 


-0.99 


894 ± 148 


111 ± 148 


8 ± 69 












0.549 


3.357 


1.91 





.66 


-0.96 


514 ± 84 


83 ± 67 


-12 ± 31 












0.550 


3.850 


2.01 





.55 


-0.88 


327 ± 59 


95 ± 54 


-52 ± 26 





.65 


± 


0. 


49 


0.550 


4.350 


2.11 





.41 


-0.83 


258 ± 48 


29 ± 62 


-13 ± 27 





.90 


± 


0. 


.43 


0.557 


4.765 


2.16 





.31 


-0.83 


222 ± 44 


-42 ± 121 


55 ± 51 


1 


.36 


± 


0. 


67 


0.601 


3.882 


1.86 





.61 


-1.26 


292 ± 57 


91 ± 89 


74 ± 50 












0.610 


4.352 


1.91 





.52 


-1.24 


221 ± 43 


75 ± 54 


-53 ± 23 












0.610 


4.850 


2.00 





.40 


-1.16 


150 ± 26 


-110 ± 48 


39 ± 20 













can only make the following conjectures as to the origin 
of this discrepancy: the Cornell results do not appear to 
have been corrected for internal virtual radiative effects 
(about 15%); their overall systematic uncertainty in ab- 
solute cross sections is 25%; their acceptance calculation 
has a model dependence which was not quantified and in 
particular the decay distribution, given in eq. © below, 
was assumed flat; finally, the estimate of average values of 
kincmatical variables (Q 2 ) and (W) may be an additional 
source of uncertainty since the corresponding bins are at 
least 5 times larger than in the present work. 



3.3 t dependence of cross sections 

Four of the 34 distributions of differential cross sections 
are illustrated in fig. El The general features of these dis- 
tributions are of a diffractive type (da/dt oc e bt ) at small 



values of — t. The values of the slope b, as determined from 
a fit of the distributions in the interval -1.5 GeV 2 < t < to, 
are between 0.5 and 2.5 GeV -2 and are compiled in ta- 
ble |3] They are also plotted as a function of the formation 
length (distance of fluctuation of the virtual photon in a 
real meson) pQ: 



V^ 2 + Q 2 + Ml - v 

in fig. El They are compatible with those obtained for 
the reaction <y*p — > p°p El- There exists only one pre- 
vious determination El °f this quantity for j*p — > tup, 
integrated over a wide kinematical range corresponding 
to 0.6 < cAt < 2.5 fm, with a value b = 6.1 ± 0.8 GeV" 2 . 
This larger value of b is consistent with the observed dis- 
crepancy between the Cornell experiment and the present 
work. For larger values of —t, the slope of the cross sec- 
tions becomes much smaller and da/dt becomes nearly 
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0.16 <x„< 0.22 0.22 < x. < 0.28 0.28 < x. < 0.34 0.34 < x„ < 0.40 




0.40 < x. < 0.46 0.46 < x„ < 0.52 0.52 < x„ < 0.58 0.58 < x„ < 0.64 
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Fig. 10. (Color online) Reduced cross sections j*p — > up as 
a function of Q 2 for differents bins in xb, in units of fib. Full 
circles: this work; open circles: Ref. |17| . The red cross and 
curves correspond to the JML model discussed in sect.0 
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Fig. 11. (Color online) Total cross sections for the reaction 
7*p — > up, as a function of Q 2 and at fixed W : this work 
in full symbols, DESY [HH (top) and Cornell [13 (bottom) in 
open symbols. Each symbol (or color) corresponds to a given 
central value of W (GeV). Note the range of integration in W 
for each data set. 
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Fig. 12. da/dt for the reaction j*p — > up, at W ~ 2.45 GeV for 
different bins in Q 2 : our data and the JML model (discussed 
in sect.QjJ with F 7T ^ 1 given by eq. 1141 (full lines) and without 
the i-dependence in this equation (dashed lines) . 



> 

o 




0.3 



0.4 0.5 0.6 

CAT (fm) 



0.7 



independent of Q 2 , except for the lowest values of W (see 
fig. I14|) . This is certainly a new finding from this experi- 
ment, which may indicate a point-like coupling of the vir- 
tual photon to the target constituents in this kinematical 
regime. This behaviour will be discussed quantitatively in 
Sect. El 



Fig. 13. Slope b of da/dt, for the reaction 7*p — > up, as a 
function of the formation length cAt, cf. eq. Q. 
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Fig. 14. da/dt at fixed values of t and W, as a function of Q 2 , 
for the reaction 7*p — > ojp : t = —0.55 (full circles), t = —1.45 
(empty circles) and t — —2.30 (squares) GeV 2 . 
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Fig. 15. (Color online) ott (open circles) and <jtl (full cir- 
cles), in units of /A>, for the reaction 7*p — > ujp as a function of 
Q 2 for different bins in xb, integrated over —2. 7 GeV 2 < t < t . 
The dashed blue and full red curves are the corresponding cal- 
culations in the JML model |29] discussed in sect.|^l 



3.4 4> dependence of cross sections 



The 34 (f> distributions have the expected <f> dependence 



da If i — - — 

— = — [a + e cos 20 utt + V 2e (l 
am 2tt \ 



e ) COS <p (TTL 



(5) 

The interference terms cttt and gtl were extracted from 
a fit of each distribution with eq. (JjJJ. The results appear 
in fig. E| and in table [3] If helicity were conserved in the 
s-channel (SCHC), these interference terms cttt and <jtl 
would vanish. It does not appear to be the case in fig. 1151 
The (j> distributions do not support the SCHC hypothesis. 



4 Analysis of u decay distribution 

In the absence of polarization in the initial state, the dis- 
tribution of the pions from lu decay is characterized by 
eq. © HID] ■ The quantities pfj are defined from a decompo- 
sition of the lu spin density matrix on a basis of 9 hermitian 
matrices. The superscript a refers to this decomposition 
and it is related to the virtual photon polarization (a = 0- 
2 for transverse photons, a — 4 for longitudinal photons, 
and a = 5-6 for interference between L and T terms). For 
example, p^ is related to the probability of the transition 
between a transverse photon and a longitudinal meson. 

All elements pfj can be expressed as bilinear combina- 
tions of helicity amplitudes which describe the 7*p — > cop 
transition |51l30| . An analysis of the W distribution can 
then be used to test whether helicity is conserved in the 
s-channel (SCHC), that is between the virtual photon and 
to. If SCHC applies, p° 00 = and p% = 1. Then eq. 



leads to a direct relation between the measured Tqq and the 
ratio R = ctl/ctt- In that case, the longitudinal and trans- 
verse cross sections may be extracted from data without 
a Rosenbluth separation. 

The matrix elements Tqq and r t ( 4 _ 1 were first extracted 
; using one-dimensional projections of the W distribution. 
I Note that r^ 4 _ 1 should be zero if SCHC applies. Integrating 



eq. © over . 
gets: 



and then respectively over <p jy or cos 9 jy , one 



W(cos M = |[d- O + (3rgS - 1) cos 2 6 N ] , (8) 



2tt 



[1 - 2r?i 1 cos2<pAr] 



(9) 



The background subtraction in the cos9n or ipjv distribu- 
tions was performed in 8 bins of the corresponding vari- 
ables, and for 8 bins in (Q 2 , xb) - The number of acceptance- 
weighted events was extracted from the corresponding Mx [epX] 
distribution, as in sect. 12.31 See figs . 1 1 61 and 1 1 71 for results, 
together with fits to eqs. © and @ respectively. 

Alternatively, the 15 matrix elements r£j- may be ex- 
pressed in terms of moments of the decay distribution 
W(cos9n, <Pn, 0) ,30. This method of expressing moments 
includes the background contribution under the u> peak 
(about 25%). It yields compatible results with the (back- 
ground subtracted) ID projection method for Tqq and rj^. 
It was used to study the t dependence of and to eval- 
uate the systematic uncertainties in their determination. 
Results for forward ^*p — > cop reaction (t! < 0.5 GeV 2 ) 
are given in fig. 1181 Systematic uncertainties originate 
from the determination of the MC acceptance. The main 
source of uncertainties was found to be the finite bin size in 
(j>. Calculations with different bin sizes (see table ^ and 
checks of higher, unphysical, moments in the event dis- 
tribution led to systematic uncertainties of 0.02 to 0.08, 
depending on the rfj matrix element. In addition, cuts in 
the event weights were varied, resulting in a systematic 
uncertainty of about 0.03 for all matrix elements. 

Finally, the rfj matrix elements were also extracted 
using an unbinned maximum likelihood method. Results 
were compatible with the first two methods. In view of the 
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W(cos6n, *Pn, 4>) = 7— ^(1 — foo) + i(3r ( ( jo — 1) cos 2 9n — v^Rer? 4 sin 29n cos tpN — r?lj sin 2 9n cos 2tpt 

47T |_ 2 2 



-e cos 2</>(ri 1 sin 2 #jv + t*qo cos 2 (9jv — \/2Rer;J sin 29n cos ipjv 
-e sm 10 sin 2#at sin (pN + Imr! _ x sin Sat sin 2<^at ) 



rj^sin 8ncos2ipn) 



+ \j2e(l + e) cos(f>(rii sin 2 9n + >"oo cos 2 9n 



V2Re 



,er 10 sin 29n cos 95 jv — r j_ x sin # jv cos 2ypjv ) 
+ y2e(l+e) sin 0(\/2Imri O sin 2#at sin (/?jv + Imrf.i sin 2 On sin 2<^at) 
where the parameters r" , hereafter referred to as matrix elements, are related to the ui spin density matrix: 



p% + eRp 



1+eR 



1 + eR 



for a = 1, 2 



1 +eJ? 



for a = 5, 6. 



(6) 
(7) 
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Fig. 16. (Color online) Distributions of acceptance-weighted 
and background-subtracted counts as a function of cos(9jv, for 
8 bins in (Q 2 , Xb)- The location and size of each graph cor- 
respond to the (Q 2 , xb) range over which the data is inte- 
grated. On all graphs, one division on the vertical axis rep- 
resents 2 x 10 4 (arbitrary units). All data are integrated in 
t (-t < 2.7 GeV 2 ). The blue curves correspond to fits with 
eq. JHJ, with the resulting Tqq and its statistical uncertainty 
indicated on each distribution. The systematic uncertainty on 
this matrix element is estimated at 0.042. 



Fig. 17. (Color online) Distributions of acceptance-weighted 
and background-subtracted counts as a function of </jjv- The 
blue curves correspond to fits with eq. ©, with the resulting 
r'l'Li and its statistical uncertainty indicated on each distri- 
bution. The systematic uncertainty on this matrix element is 
estimated at 0.042. See also legend of fig. 1161 
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4> dependence of the acceptance (see fig. El, this method 
was used for checking the validity of the rfj determination 
when restricting the <p range taken into consideration in 
the fit. 

These studies lead to the conclusion that SCHC does 
not hold for the reaction 7*p — > cup, not only when con- 
sidering the whole t range (fig. I17JI. but also, though in a 
lesser extent, in the forward direction (fig. lTH)) . For SCHC, 
all matrix elements become zero, except five: 
Im rf_ l5 Re r 
they satisfy 0. 
The quantity 



r 00> ' 1-1) 

f , Im rf and these are not all independent; 

1 Tth anrl I -) i-i(ji5 



= — Imrf_ 1 and Re? 
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= -Imr? . 
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E 



r 

Ar 



(r 1 ^ + lmrl_ x ) 2 
(zAr 1 .,) 2 + (Almrl^y 



where the sum is carried over the ten matrix elements 
which would be zero if SCHC applies, may be used as a 
measure of SCHC violation. Including in the denomina- 
tors Ar the systematic uncertainties added in quadrature 
to the statistical uncertainties, the 7 x 2 values (excluding 
the distributions at the lowest xb bin where SCHC viola- 
tion is the most manifest in fig. 118(1 range from 2.3 to 7.7 
when including all data, and drop only to 1.7 to 5.1, in 
spite of doubled statistical uncertainties, when consider- 
ing only the forward production (t' < 0.5 GeV 2 ). Further- 
more, when examining the relation between these matrix 
elements and helicity-flip amplitudes, it does not appear 
possible to ascribe the SCHC violation to a small subset 
of these amplitudes. It is therefore not justified to cal- 
culate R from eq. Q and separate the longitudinal and 
transverse cross sections from this data. 
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Fig. 18. (Color online) rfj extracted with the method of mo- 
ments for 8 bins in (Q 2 , xb) and for t 1 < 0.5 GeV 2 . The loca- 
tion and size of each graph correspond to the (Q 2 , xb) range 
over which the data is integrated, but the scale is the same 
on all graphs. The abscissa on each graph corresponds to the 
following list of matrix elements: roo, Re r?o, r?*i, r^o, rji, 
Re rlo, r{-i, Im r 2 , Im r 2 ^, rgo, rfi, Re rfo, rf_i, Im rf , 
Im rf_!. The red filled symbols indicate those matrix elements 
which are zero if SCHC applies. The 16th entry (blue empty 
circle, in some cases off scale) is the combination of rfj given 
by eq. Ulljl . Error bars include systematic uncertainties added 
in quadrature. 



When one retains only those amplitudes which corre- 
spond to a natural parity exchange in the i-channel, then 
the following relation should hold : 



1 



„04 



2r» 
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2r\ x - 2ri_! = 



(11) 



This particular combination is plotted as the 16th point 
on each of the graphs of fig. El The fact that it is not zero 
points to the importance of the unnatural parity (presum- 
ably pion) exchange. 

It is also possible to estimate qualitatively the role of 
pion exchange through the U/N asymmetry of the trans- 
verse cross section, where U and N refer to unnatural and 
natural parity exchange contributions 5 : 



P = 



_ of 



= (l+eR)(2r 1 1 _ 1 -r 1 m ) 



(12) 



Our results yield r\_ x < and Tq 
kinematical range, and thus: 

P<-(2|ri_ x | ' ; 



> over the whole 



oo 



I) 



(13) 



Hence P is large and negative, which means that most 
of the transverse cross section is due to unnatural parity 
exchange. 



5 Comparison with a Regge model 

Regge phenomenology was applied with success to the 
photoproduction of vector mesons in our energy range 



Table 4. Meson and Pomeron (or two-gluon) exchanges con- 
sidered in the JML model for vector meson production. 



Produced 


Exchanged 


vector meson 


Regge trajectories 


P 


e, h, P/2g 


Id 


7T°, h, P/2g 




P/2g 



and at higher energies |4llST] . Laget and co-workers showed 
that the introduction of saturating Regge trajectories pro- 
vides an excellent simultaneous description of the high 
—t behaviour of the jp — > pp, w, 4> cross sections, given 
an appropriate choice of the relevant coupling constants. 
The i-channel exchanges considered in this JML model are 
indicated in table 01 Saturating trajectories have a close 
phenomenological connection to the quark- antiquark in- 
teraction which governs the mesonic structure |221. They 
provide an effective way to implement gluon exchange be- 
tween the quarks forming the exchanged meson. 

This model was extended to the case of electroproduc- 
tion j2H]. The Q 2 dependence of the f2 and P exchange 
is built in the model. In the case of lu production, the 
only additional free parameters come from the electro- 
magnetic form factor which accounts for the finite size 
of the vertex between the virtual photon, the exchanged 
7T° trajectory and the u) meson. This form factor could 
be chosen as the usual parameterization of the pion elec- 
tromagnetic form factor: F W7r7 = F n = (1 + Q 2 / A 2 .)^ 1 , 
with A 2 = 0.462 GeV 2 . As described so far, the model 
fails to account for the observed t dependence (see dashed 
lines in fig. H2f> . From the observation that the differential 
cross section becomes nearly Q 2 -independent at high —t, 
an adhoc modification of the form factor 



1 



( l±2i 



(*) 



(°> 



(14) 



was proposed [2H|- The saturating n° Regge trajectory 
obeys the relation lim t ^_oo a^(t) = —1, so that the form 
factor becomes flat at high —t. Thus, eq. (|14f> associates 
the point-like coupling of the virtual photon with the satu- 
ration of the 7T° Regge trajectory which accounts for hard 
scattering in this kinematical limit [2S| - Note that this 
modification of the form factor does not violate gauge 
invariance, which holds separately for each contribution 
from Table 01 and, in the case of ir° exchange, from the 
spin and momentum structure of the wk^/ vertex. 

The t dependence of the differential cross sections is 
then well described (solid lines in fig. 1121) . The Q 2 de- 
pendence of the cross sections is illustrated in fig. ^3 At 
high xb, which corresponds to the lowest values of W, 
s-channel resonance contributions are not taken into ac- 
count in the model and may explain the observed disagree- 
ment. Finally the interference terms cttt and otl agree 
in sign and trend, but not in magnitude, with our results 

(fig. El. 
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So, within this model, tt° exchange, or rather the ex- 
change of the associated saturating Regge trajectory, con- 
tinues to dominate the cross section at high Q 2 and the 
cross section is mostly transverse. This is consistent with 
our observations of the dominance of unnatural parity ex- 
change in the i-channcl in the previous section. 



6 Relevance of the handbag diagram 

Let us recall that the handbag diagram of fig. His expected 
to be the leading one in the Bjorken regime. In this picture, 
the transition 7^ — > ujl would dominate the process. This 
is clearly antinomic to the findings in sect. where our 
results are interpreted as dominated by the tt° exchange, 
which is mostly due to transverse photons. In addition, ir° 
exchange is of a pseudo-scalar nature, while the H and E 
GPD which enter the handbag diagram amplitude are of 
a vector nature. 

Independent of the model interpretation presented in 
sect.0 our results point to the non-conservation of helicity 
in the s-channel ffigs. 1151 1171 and I18fl . meaning that the 
handbag diagram does not dominate the process, even for 
small values of —t and Q 2 as large as 4.5 GeV 2 . 

As a consequence, ol could not be extracted from our 
data for a direct comparison with models based on the 
GPD formalism. It is however instructive to consider here 
the predictions of a GPD-based model [5055) . denoted 
hereafter VGG. This is a twist-2, leading order calcula- 
tion, where the GPD are parameterized in terms of dou- 
ble distributions (DD) and include the so-called D-term 

(see Ref. jSS] for definitions): H,E ~ DD(x, ^e 6 ^ 3 )*/ 2 , 
where b is taken from the data (see sect. 13.31 and tablc|HJ). 
An effective way of incorporating some of the higher twist 
effects is to introduce a "frozen" strong coupling constant 
as = 0.56. This model is described in some more details 
in Ref. fT^j and is applied here to the specific case of ui 
production. The model calculations (VGG and JML) of 
£<tl are plotted in fig.f^l The sharp drop of the curves at 
high Q 2 is due to the decrease of e, at our given beam en- 
ergy, as Q 2 reaches its kinematical limit. When compared 
to our results, eat is calculated to be only 1/6 to 1/4 of 
the measured cross sections, thus explaining the difficulty 
in extracting this contribution. 

The u> channel thus appears to be a challenging reac- 
tion channel to study the applicability of the GPD formal- 
ism. This is attributed to the i-channel ir° exchange, which 
remains significant even at high values of Q 2 . In contradis- 
tinction, the 7T° exchange is negligible in the case of the p 
production channel, where SCHC was found to hold, and 
<jl could be extracted and compared successfully to GPD 
models |TT1T2| . 



7 Summary 

An extensive set of data on exclusive u> electroproduction 
has been presented, for Q 2 from 1.6 to 5.1 GeV 2 and W 
from 1.8 to 2.8 GeV (x B from 0.16 to 0.64). Total and 




2 3 
Q 2 (GeV 2 ) 



Fig. 19. (Color online) Total cross sections for the reaction 
7*p -> cop, for (W) = 2.1 (top) and 2.8 (bottom) GeV : this 
work (full red circles), DESY data (empty diamonds), Cornell 
data (empty circles), and JML model (dotted curves). The lon- 
gitudinal contribution, e(E,Q 2 )aL, is calculated according to 
the JML (solid lines) and VGG (dashed blue lines) models. 



differential cross sections for the reaction 7*p — > up were 
extracted, as well as matrix elements linked to the u> spin 
density matrix. 

The t differential cross sections are surprisingly large 
for high values of —t (up to 2.7 GeV 2 ). This feature can 
be accounted for in a Regge-based model (JML), provided 
a t dependence is assumed for the lottj vertex form factor, 
with a prescription inspired from saturating Regge trajec- 
tories. It appears that the virtual photon is more likely to 
couple to a point-like object as —t increases. 

The analysis of the <j> differential cross sections and of 
the uj decay matrix elements indicate that the s-channel 
helicity is not conserved in this process. As a first conse- 
quence, the longitudinal and transverse contributions to 
the cross sections could not be separated. Furthermore, 
the values of some decay matrix elements point to the im- 
portance of unnatural parity exchange in the i-channel, 
such as 7T° exchange. This behaviour had been previously 
established in the case of u> photoproduction, but not for 
the large photon virtuality obtained in this experiment. 
The results on these observables also support the JML 
model, where the exchange of the saturating Regge tra- 
jectory associated with the n° is mostly transverse and 
dominates the process. 



14 



Finally, the experiment demonstrated that exclusive 
vector meson electroproduction can be measured with high 
statistics in a wide kincmatical range. The limitations at 
high Q 2 were not due to the available luminosity of the 
CEBAF accelerator or to the characteristics of the CLAS 
spectrometer, but to the present beam energy. With the 
planned upgrade of the beam energy up to 12 GeV |M4| . 
such reactions will be measured to still higher values of Q 2 . 
In the specific case of the to meson, as was shown in this 
paper, this will be a necessary condition for the extrac- 
tion of a longitudinal contribution of the handbag type, 
related at low values of —t to generalized parton distribu- 
tions. More generally, this experiment opens a window on 
the high Q 2 and high —t behaviour of exclusive reactions, 
which needs further exploration. 
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